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The effects of a radiation field (RF) on the unstable modes developed in relativistic electron beam-plasma 
rvj , interaction are investigated assuming that cjq > ijjp, where ujq is the frequency of the RF and Up is the plasma 

frequency. These unstable modes are parametrically coupled to each other due to the RF and are a mix between 
f^ ' two-stream and parametric instabilities. The dispersion equations are derived by the linearization of the kinetic 

^Sl . equations for a beam-plasma system as well as the Maxwell equations. In order to highlight the effect of the 

> . ' radiation field we present a comparison of our analytical and numerical results obtained for nonzero RF with 

C^ ' those for vanishing RF. Assuming that the drift velocity Ub of the beam is parallel to the wave vector k of the 

excitations two particular transversal and parallel configurations of the polarization vector Eq of the RF with 
respect to k are considered in detail. It is shown that in both geometries resonant and nonresonant couplings 
1/^ ' between different modes are possible. The largest growth rates are expected at the transversal configuration 

^SJ , when Eo is perpendicular to k. In this case it is demonstrated that in general the spectrum of the unstable modes 

in u)~k plane is split into two distinct domains with long and short wavelengths, where the unstable modes are 
mainly sensitive to the beam or the RF parameters, respectively. In parallel configuration, Eo || k, and at short 
wavelengths the growth rates of the unstable modes are sensitive to both beam and RF parameters remaining 
insensitive to the RF at long wavelengths. 

PACS numbers: 52.40.Mj, 52.40.Db, 52.35.Mw, 52.35.Qz 
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O . I. INTRODUCTION 

• 1—1 
^ . 

^^' The interaction of a relativistic electron beam (REB) with a plasma is a subject of relevance for many fields of physics rangin: 

from inertial fusion to astrophysics |l-5]. This interaction is also relevant, among others, for fast ignition scenario (FIS) 

where the precompressed deuterium-tritium (DT) core of a fusion target is to be ignited by a laser-generated relativistic electron 

beam. The REB quickly prompts a return current so that one eventually has to deal with a typical two-stream configuration 

which is subjected to various electromagnetic instabilities. Much effort have been devoted in the last years to investigate these 

instabilities 18- 211 . whether it be the two-stream 11221 12311 . the filamentation |24] or Weibel 1251 instabilities. These instabilities 

are usually analyzed by linearizing the relativistic Vlasov or fluid and Maxwell equations. Then, the response of the linearized 

equation to a perturbation is investigated and one eventually finds some unstable self-excited modes. At this stage, the orientation 

(V> of the wave vector k of the excitations plays an important role. Choosing a wave vector parallel to the beam velocity Uf, yields 

• the two-stream unstable modes which are of electrostatic nature and therefore characterized by wave and electric field vectors 

\^ both parallel to the beam propagation direction. On the other hand, choosing a wave vector normal to the beam yields the purely 

*vj transverse filamentation unstable modes. They are mostly electromagnetic, purely growing, and develop preferentially in the 

,_l ■ plane normal to the beam. The exploration of the much less investigated intermediate orientations has brought a very important 

i;^ I result b y th at showing that the strongest instability suffered by the system is eventually to be found for an oblique wave vector 

III4I - I20I l26l [2711 ■ This most unstable mode is a mixture of the two-stream and the filamentation instabilities but it is not damped 

as the last two ones. For example, the maximum two-stream growth rate is reduced by a factor of 7^^ in the relativistic regime, 

^ while the most unstable mode only decreases by a factor of 7^^/'^, where 7 is the beam relativistic factor. The filamentation 

growth rate varies as 7^^/^ and may be reduced, even canceled, by a transverse beam temperature 1I12I1 whereas the most unstable 

mode is quite insensitive to temperatures as long as they are nonrelativistic [15]. 

In general beam-plasma instabilities have been studied in detail for many physical situations including the interaction of the 
cold, warm, inhomogeneous or anisotropic electron or ion beams with cold, warm, magnetized, unmagnetized, inhomogeneous 
or anisotropic plasmas, see, e.g., Refs. [^-'21'] and references therein (see also Refs. |28, 29] for detailed bibliography). In this 
paper, we present a study of the effects of a radiation field (RF) on the interaction of REB with a plasma. More specifically our 
objective is to study the beam-plasma instabilities in a laser irradiated plasma which, to our knowledge, has not been discussed 
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in the literature. The principal motivations of the present paper are the researches on the topic of the FIS for inertial confinement 
fusion li6, J] which involves the interaction of a laser-generated REB with a hot plasma. Although the concept of the FIS implies 
an overdense plasma and the propagation of a relativistic electron beam from the border of a pre-compressed target to the dense 
core occurs without crossing the laser beam, the target plasma is assumed to be parametrically excited by the RF through high 
harmonics generation. In addition, a promising inertial confinement fusion scheme has been recently proposed 130-321 . in 
which the plasma target is irradiated simultaneously by intense laser and ion beams. In both situations the presence of the RF 
can dramatically change the main features of the standard (i.e. when the laser is off) beam-plasma instabilities. 

Previously the interaction of charged particles with plasma in the presence of the RF has been a subject of great activity in the 
contexts of the stopping power 1133143611 (see also references therein) and other processes of interest for applications in optics, 
solid-state, and fusion research. In particular, the analytical calculations Il33l - [35ll as well as the particle-in-cell (PIC) numerical 
simulations 1 36] show that the propagation and the subsequent stopping of the charged particles would be essentially affected by 
the parametric excitations of the plasma target by means of laser irradiation. It is well known that in general in the absence of the 
charged particle beams the laser irradiated plasma is subjected to the parametric instabilities 11371 13811 . Therefore in the present 
context it is expected that the above mentioned beam-plasma unstable modes developed in a plasma irradiated simultaneously 
by laser and electron beams are parametrically coupled to each other and are not, in principle distinguishable. 

The plan of the paper is as follows. In Sec. [Ill we outline the kinetic formulation for the interaction of the REB with a laser 
irradiated plasma. The full electromagnetic response of the plasma is considered. The general linear dispersion relations are 
derived in Sec.|IIIl which are studied for two particular cases in Sees. |IV| and [Vl with transversal and parallel configuration of 
the polarization vector of the RF with respect to k, respectively. Furthermore, assuming nonrelativistic laser intensities only the 
lowest (zero, first and second) harmonics of the electromagnetic fields are considered. The obtained dispersion equations are 
investigated numerically in Sec. |VT] The results are summarized in Sec. IVIII which also includes discussion and outlook. In 
Appendix lAl we consider the standard (in the absence of the RF) stable and unstable modes of the beam-plasma system in a 
cold-fluid approximation. The asymptotic behavior of the frequencies and the growth rates of these modes at large and small 
k are considered. In Appendix iBl we provide some technical details for an evaluation of the sum containing Bessel functions. 
An equation describing the evolution of the amplitude of the parametrically excited plasma waves and involving all excited 
harmonics is derived and discussed in Appendix |C] 

II. THEORETICAL BACKGROUND 

In this section, we consider the main aspects of the interaction of the relativistic electron beam (REB) with a homogeneous, 
collisionless, and unmagnetized plasma in the presence of high-frequency radiation field (RF), Eo(i) = Eosin(a;ot). Here 
Eo and cuo are the amplitude and the frequency of the RF. The problem is formulated using the perturbative treatment, and 
includes the effects of the RF in a self-consistent way. The RF is treated in the long-wavelength limit, and the plasma particles 
(electrons and ions) are considered nonrelativistic. These are good approximations provided that (1) the wavelength of the RF 
(Aq = 27Tc/uJo) is much larger than the Debye screening length Ad = Vth/<^p with wth the thermal velocity of the electrons and 
ujp the plasma frequency, and (2) the "quiver velocity"of the electrons in the RF (ve — eEo/moJo) is much smaller than the 
speed of light c. These conditions can be alternatively written as (1) ujo/ojp <C 27rc/t;th, (2) II ^ ^"■e"^c'^('^o/^p)^' where 
II ~ cEq/Stt is the RF intensity. As an estimate in the case of dense plasma, with electron density Ue — 10^^ cm^"^, we get 
ingTTic^ ~ 1.2 X 10^° W/cm^. Thus the limits (1) and (2) are well above the values obtained with currently available high-power 
RF sources, and so the approximations are well justified. Furthermore we consider an underdense plasma with ojq > tUp. 

The dynamics of the beam-plasma system is governed by the relativistic and nonrelativistic Vlasov kinetic equations for 
the distribution functions of the REB fi,{r,v,t) and the plasma fa{r,v,t) (where a — e,i indicates the plasma species), 
respectively, as well as by the Maxwell equations for the electromagnetic fields. Thus, 

^^+v.^-+e„,,|Eo(t)+E + -[vxB]|.^-=0, (1) 

where Ca and et are the charges of the plasma and beam particles, respectively, E and B are the self-consistent electromagnetic 
fields which are determined by the Maxwell equations. 
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yUa + CbTlfc I , V • B = 0. (3) 

Here na{v,t) and nb{r,t) are the densities for the plasma species a and for the relativistic beam, respectively, JQ(r,i) and 



jb(r. t) are the corresponding currents induced in plasma and beam, respectively, 

na,b (r, *) = / fa,b (r, P, t) dp, ja,b (r, t) = e^^b / ^fa,b (r, P, t) dp. (4) 

As mentioned above we consider a nonrelativistic plasma and in Eq. ([T]i for the distribution function /„ the momentum 
connects linearly to the particle velocity, p — TOqV, where nia is the mass of the plasma species a. Equation ([T]i for the 
distribution function /;, is relativistic and p = m-b^fv in this case, where 7 = (1 — w^/c^)~^/^ and rrib are the relativistic factor 
and the rest mass of the beam particles. Moreover, we consider an ultrarelativistic electron beam with 7^ — (1 — u^/c^)^^/^ ^ 1, 
where Ub is the beam drift velocity, and therefore the influence of the RF Eo(i) on the beam is ignored in the kinetic equation 
([T]l for the distribution function /b(r, p, t). 

For sufficiently small perturbations, we assume fa.b = foa.b + fia.b (with fia.b ^ foa.b), where /qq and fob are the 
equilibrium distribution functions of the plasma species a and the beam in an unperturbed state, respectively. The solution of the 
linearized kinetic equation ([T]i for the relativistic beam, when the RF Eo(i) in Eq. ([T]i is neglected, is well known. This standard 
calculation is explained at length in a number of plasma physics textbooks, [38-41], and we just here mention the final result. In 
terms of the Fourier transformed quantities the solution of the kinetic equation reads 



fib (k, uj, p) = iebEj (k, io) 
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Here E(k, uj) and /i6(k, uj, p) are the Fourier transforms of the electric field and the beam distribution function, respectively, 
with respect to variables r and t, Sij is the unit tensor of rank three. Note that the positive infinitesimal -t-zO in Eq. Q guarantees 
the causality of the response. 

The perturbations of the densities and the currents induced in the plasma and the beam are determined from Eq. (|4). The 
Fourier transforms of these quantities are then given by 



ni„,6 (k, uj) ^ fia,b (k, UJ, p) dp, ji^^b (k, uj) = ea,b / ^fia.b (k, uj, p) 



dp. 



Substituting the distribution function (|5]l into these quantities and integrating by parts with the help of the relation dvi/dpj 
{l/mb'j){5ij — ViVj/c^) (where 7^(p) = 1 + p'^/mlc^) for the induced current and the density of the beam we obtain 



jib,i (k, uj) = ab.ij (k, uj) Ej (k, uj) , 
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where ab^ij (k, ui) is the conductivity tensor of the relativistic charged particle beam. 
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Consider now the solution of the kinetic equation ([T]i for the plasma electrons and ions in the presence of the high-frequency 
RF. In an unperturbed state (i.e. neglecting the self-consistent electromagnetic fields E and B in Eq. ^ and assuming the 
homogeneous initial state) the distribution function satisfies the equation 



dfi 



Oq 



CaEo sin {ujot) 



dfoo 







which yields the equilibrium distribution function for the plasma species a, 

foa (P, t) = F„ (p + rUaVEa COS {ujQt)) . 



(10) 



(11) 



Here Fa{p) is an arbitrary function. Below we will assume that this function is isotropic in momentum space, vsa — 
eaEio/iriaUJo, a^ = eaEo/niaUiQ are quiver velocity and the oscillation amplitude of the plasma species a, respectively, driven 
by the RF. 

Next we consider the linearized kinetic equation for the plasma species a. Introducing the Fourier transforms /ia(k, p, t), 
E(k, t) and B(k, t) with respect to the coordinate r the linearized kinetic equation reads 

(12) 
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dp ■ 



In order to solve Eq. (fT2] l it is convenient to introduce a new unknown function ^^ via relation 

/ia(k,p,t) = e<"^'"("°*)vl/„(k,P,i), 



(13) 



where Cq = (k • a^^) and P = p + niaVEa cos{ujf)t). Substituting this relation into Eq. (fT2] i it is easy to see that the ob- 
tained equation for the unknown function 'ia constitutes an equation with periodic coefficients where the role of the dynamic 
momentum is now played by the quantity P. Therefore we introduce the decomposition 
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-OO ^^ 
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(14) 



where (5(k, p, t) represents one of the quantities ^^(k, p,i), E(k, t) and B(k, i) and Q'^")(k, w,p) is the corresponding am- 
plitude of the nth harmonic. From Maxwell equation we express the magnetic field through the electric field. In terms of the 
amplitudes of the nth harmonics (see definition given by Eq. (fl4] i) this relation is given by 
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Also in the kinetic equation derived for \l/a ' taking into account the isotropy of the equihbrium distribution function Fa of the 
plasma particles we neglect the term 



dFa 

[hvj - % (k • v)] • -^ = 



(16) 



and using the Fourier series representation of the exponential function ll42ll (see also Eq. ( IC2l l) for the amplitude of the nth 
harmonic of the distribution function we finally obtain 



vl/i")(k,w,p) = - 

X ^ Jn-liCa) 



dFa (p) 
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(17) 
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Here Jn is the Bessel function of the nth order. Throughout this paper the upper indices given in the parentheses indicate the 
harmonic number while the lower indices determine the components of the vectors and tensors. 

The amplitudes of the harmonics of the induced current and the charge density in a plasma are obtained from Eqs. (|6}, ( fTJt . 
(O and (fTTb . Straightforward calculations yield 
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for the Fourier transforms of the sth harmonics of the current and the charge density, respectively, with 

k u! + nojQ Qa 
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(20) 



Here we have introduced the conductivity tensor a)- (k, uj 



_, oo 



(o; + nuj^Y 

a n— — oo 



^n—i t,S>aj ^n — s \SctJ 



(21) 



th r, th'p 



'jr 



Sea± (n) + xf "^ (k, lo) xi''^^ (k, oj) 6s ^n (n) 



and the abbreviations Sea\\:±{'n') = Sea\\:±{k, uj + nujo), where 5ea\\ {k, uj) and 6ea±{k, tu) are the partial contributions of the 
plasma species a to the longitudinal and transversal (with respect to the wave vector k) dielectric functions (see, e.g., |39]), 
respectively. 
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Note that since the equilibrium distribution function F^ is isotropic the partial dielectric functions Ssa || ; j^ (fc, w) are also isotropic, 
i.e. they do not depend on the direction of the wave vector k. The obtained expressions (|T8]l-(|23]l with Eqs. (jT)-® as well as 
the Maxwell equations (|2|i and (O written in the Fourier space completely determine the electromagnetic response in the beam- 
plasma system in the presence of the RF. Using this system of equations the general dispersion equations are derived in the next 
section. 

We would like to close this section with the following two remarks. First, the distribution function /i6(k, w, p) of the beam 
particles given by Eq. (|5]l as well as the induced current (|7) and the charge density (|8]l are determined by the Fourier transform 
E(k, oi) of the electric field. In contrast to this case the distribution function (Eq. ( fT3] l with Eqs. ( fT4l i and (fTTji), the induced 
current (Eq. (fT¥t ) and the density (Eq. ([19)) of the plasma are determined by the Fourier transform of the amplitude of the 
harmonics of the electric field, E'^"^(k, w). From Eq. ( fT4] i it is straightforward to deduce the connection between the Fourier 
transforms E(k.cj) and E'^"'(k, w). Changing the integration variable in each term of summation in Eq. ( fT4b according to 
UJ + nujQ — >■ a; we obtain 



E(k,o.) 



E 



E(") (k, cj - nwo) . 



(24) 



Thus, E(k, w) is the sum of all harmonics E^"' (k, cj) with shifted frequencies w ± nwo- Second, assuming an ultrarelativistic 
beam for derivation of the distribution function (|5} we have neglected the RP in the kinetic equation ([T]) for /{,. And as a 
consequence the perturbation of the beam distribution function is determined by E(k.ci;). Although the RF is not directly 
involved in the kinetic equation ^ for /b, it affects this distribution function via the self-consistent electric field E(k,(u;) 
containing all harmonics produced by the RF (see Eq. (|24]|). 



III. DISPERSION EQUATION 



In this section using the expressions derived for the induced currents in the beam and plasma we consider the dispersion 
equation of the waves excited in a plasma by the relativistic beam of charged particles. For this purpose we employ the Maxwell 
equations ^. Introducing Fourier transforms of the electric field and the currents according to Eq. ([14) and excluding the 
magnetic field from these equations by means of Eq. ( [TST i from Eqs. ([2]l, ([2ll-([9|) and ([T8]l-(i2ni for the components of the 
amplitude of the nth harmonic of the electric field we obtain 
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Here the conductivity tensors of the plasma cry" (k, uj) and the beam ab^rj (k, a;) are determined by Eqs. (1211 1 and ([9), respec- 
tively. It is seen that in the right-hand side of Eq. ( |25l l the beam current vanishes for any nonzero harmonic number, n ^ 0. 

Before starting the systematic investigation of the general dispersion equation for the beam-plasma system in the presence of 
the RF it is constructive consider briefly two limiting cases of Eq. (i25] l. First, at the vanishing RF (i.e. at (a ~> 0) from Eq. (1211 1 



{nl). 



it is straightforward to calculate the conductivity tensor ct, (k, a;) of the plasma which reads in this limit 



ry ' iKuj) = dni ^^. [e^-j (n) - d^j\ 



SnlCTij (n) . 



(26) 



Here we have introduced (as above) the abbreviations (Jij{n) = (7ij(k,a; + nujo), eij{n) = £ij(k,aj + nuj^), £||;_L(n) = 
e||.^(fc, UJ + nujf)), and CTij(k, uj) and 



e^j (k,w) = -y-eii {k,uj) + Sij - -—f- e^ {k,uj) 
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(27) 



are the conductivity and the dielectric tensors of an isotropic plasma, respectively, with longitudinal (e\i {k, uj)) and transversal 
(e^{k, ui)) dielectric functions (see, e.g., Refs. 11381 [39ll ). 

£l±{k,uj) = l + ^6eai±ik,uj). (28) 

a 

Substituting Eq. ( l26l ) into Eq. dZSl) we obtain that £^ (k, uj) — SnoEj (k, u) and Eq. (|24] | is fulfilled automatically. Thus, after 
some simplifications we arrive at 

jfc^jy - hk, - ^ [£„■ {k,oj) + Set.j {k,uj)]\ Ej (k,c^) = 0. (29) 

Here £ij (k, oj) + Seb.ij (k, w) is the total dielectric tensor of the beam-plasma system, where Seb.ij (k, i^) = (47ri/w)crfc,y (k, oj) 
is the partial contribution of the beam to the total dielectric tensor of the system. Equating the determinant of the system of 
linear equations (|29] l to zero yields the general dispersion equation for the beam-plasma system. In the last years this equation 
has been studied in detail for arbitrary orientation of the electron beam propagation direction with respect to the wave vector k, 
see, e.g., Refs. lil0l - l2lll . In particular, assuming cold, homogeneous and monochromatic charged particle beam with unperturbed 
distribution function /ofc(p) = nfc(5(p — pb), where rn, is the beam density, from Eq. (|9]l one obtains 



5eb,ij (k, uj) ■ 



JbUI^ 
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(30) 



Here pf, = TUb^/bUb, 7fc = (1 — u^/c^)~^/^ and o;^ — innbel/mb are the relativistic factor and the plasma frequency of 
the beam, respectively. At this stage it is convenient to represent the vectors (particularly the electric field) in the form of an 
expansion in the components parallel, Ay = (k • A)/fc, and perpendicular, A^ = A — (k/fc)A||, to the wave vector k. In 
particular, choosing a wave vector k parallel to the beam velocity Ub yields the two-stream (TS) unstable modes which are of 
electrostatic nature with E^^(k,uj) ^ and E j^ (k, w) — 0. Introducing the longitudinal dielectric function of the beam by means 
of the relation, 

Ssbw (k,^) = ^Ssb.j (k,^) = ~——^ ^2' (31) 

the dispersion equation in this case for a beam-plasma system then reads 

X>|| (fc, W) = £|| (fc, LU) + fcfcll (fc, LO) = 0. (32) 

On the other hand, choosing k normal to the beam velocity ut, yields the purely transverse (electromagnetic) filamentation un- 
stable modes. It should be emphasized that we have considered above an infinite beam of charged particles and as a consequence 
the plasma return current is not involved in Eqs. (|3CTt - (l32i in self-consistent manner This is not, however, a strong limitation 
of the present treatment. For instance, the drift velocity Ue of the plasma return flow can be deduced from the beam current 
neutralization condition, n^Ue ~ —UbUb. Then within a cold-fluid model the return current is included by adding in Eqs. ( l30l ) 
and ( [3TT i the similar terms but with a flow velocity Ug and plasma density Ug llTJl - Uql . 

Second, in the case of the absence of the external beam the last term in the right-hand side of Eq. ( |25t vanishes. Then the 
remaining infinite system of equations for the electric field harmonics represents the electromagnetic response of the plasma to 
the RF. In general the longitudinal and transversal components of the electric field are coupled parametrically and the excitations 
are a mixture of both types of modes. Previously the parametrically unstable modes have been studied in detail both for the 
electrostatic ll43ll (see also Refs. |37, 38]) and electromagnetic |44, 45] excitations. The purely electrostatic excitations with 
E^^ = are possible when the polarization vector of the laser radiation is parallel to the wave vector k (Eo± — 0). In this case 
the plasma electrons and ions are driven by the laser field only in the direction of k. 

To illustrate the problem of the charged particles beam-plasma instabilities developed in a laser irradiated plasma, we consider 
below two examples when the polarization vector Eq of the RF is perpendicular (Sec. IIV) or parallel (Sec. |V} to the wave vector 
k. We consider an infinite and cold beam of charged particles of velocity Ub aligned with the direction of k and uniform density 
Ub passing through a homogeneous electron plasma with density of electrons n^. Therefore the partial contribution of the beam 
to the total dielectric function of the beam-plasma system is given by Eq. ( [30l l. For simplicity we will use throughout the 
notation Ub\\ = Ub for the beam velocity. In the case the RF is off the chosen geometry corresponds to the excitations of the TS 
unstable modes provided that the return plasma current is included in the dispersion relations. Nevertheless, in the present study 
neglecting the return current we will adopt the terminology "two-stream instability". This should not be confusing as long as 
the velocity of the beam is parallel to k. 



IV. TRANSVERSAL POLARIZATION OF THE RE (Eq ± k) 

In this section we consider Eq. dZST l for the harmonics of the electric field when the polarization vector Eq of the laser field is 
perpendicular to k((kEo) =0 and Cq, = 0). Then from Eq. (ISTT i it is seen that the nonvanishing components of the conductivity 
tensor are a^-' ', a^.-' , cr^ ' ' ^ while a^-' ' = at p ^ 3. Using this fact Eq. ( |25l l for the electric field harmonics is 
represented as 
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where we have introduced the following notations E^") (i) = E^") (k, w - Il^q), E(^) ee E(k, w - £wo), S^j (t) = T,rj (k, 
£uja), D^^{i) = D^^{k,uj + £ujo), R%{i) = R%{'k,Lo + lua), feb,.jW = 5£6^,j(k, w + ^wq) and 
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In the following we consider throughout an electron plasma neglecting the dynamics of plasma ions. Thus we restrict ourself 
by the frequency domain well above the ionic frequencies. We introduce the oscillation amplitude and the quiver velocity of the 
electrons via relations a^ = — ae, wec = —ve^, ve = eEo/mujQ, a = eEo/mujQ (— e is the electron charge), and e ~ Eq/Eq. 
So, the quantities a and ve are positive by definition. 

Next, for exclusion of harmonics E- (0) in Eq. (l33T l the frequency uj is replaced by a; — iojQ and using Eq. ( l24b we perform 
summation over £. This yields an equation for the amplitude Ej (0) of the electric field, 

,2 



k^Sr 



rvp hj rt 



^S,,(0) 



E,{Q) 



= D+{0)E,{-1) + D;^{0)E,{1) + R+{0)E,{~2) 
+R;A0)E,i2) + ^5eb^rA^)EM- 



(37) 



The resulting equation represents an infinite system of linear equations for the quantities Ej{ztp) (with p = 0,1, 2, ...). The 
(infinite) determinant of this system determines the dispersion equation for the beam-plasma system in the presence of the RF. 
It follows from Eq. ( |37] | that for the perturbations of which electric vector E(0) is polarized perpendicular to the plane of the 
vectors k and Eg, there is no instability. The dispersion equation for these modes is given by 'D± {k,uj) = 0, where 



V^ik,uj) = l 



1 



fc2c 



^^u;^e^ik,uj) 
lb 



(38) 



This is simply the dispersion equation for the ordinary transverse electromagnetic modes propagating in an isotropic plasma 
0381 13911 in the absence of the RF but modified due to the presence of the beam. The second term in the right-hand side of 
Eq. ( [38l l with minus sign is the partial contribution of the transverse dielectric function of a cold beam to the total transverse 
dielectric function of the beam-plasma system. It is noteworthy that this contribution depends on 7^^ while the longitudinal 
contribution ( [3T] i decays as 7^^ with the beam relativistic factor 

To reveal the instability we consider therefore the case of a polarization wherein the electric vector E(0) of the perturbations 
lies in the plane containing the vectors k and Eq. Introducing the components of the electric field parallel (_E||) and perpendicular 
(E^) to the wave vector k, for these modes from Eq. d37T i we obtain 



vmEM) 



kvE 



[^(l)+V(-l)]fc||(0), 



(39) 
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^^(0)-74[fe||(l) + <5£||(-l)] ^(0) 



4c2 



vt 



4c2 

ve 



[i^(-2)fc||(l) + ^(2)fc||(-l)] (40) 



+ :^[^ll(-l)'^^ll(l) + ^ll(l)'5^ll(-l)]' 



where X'ii(^) = Vi\{k,uj + £ujQ),Vj_{e) = Vj_{k,Lu + ecuo), andip{k,uj) = (e • E^(k, w))/a; with ^(^) = V(k,a; - ^wq). Let 
us recall that the function 2?|| (0) given by Eq. ( [32] i is the total longitudinal dielectric function of the beam-plasma system in the 
absence of the RF. In this case the transverse and longitudinal modes are independent with the dispersion relations V^ (0) = 
and I?|| (0) = 0, respectively. However, in the presence of the laser radiation these modes are parametric ally coupled according 
to Eqs. (3% and (l40l i. The longitudinal electric fields in Eq. ( l40l l can be excluded inserting the values _E||(— 1) and £^||(1) 
calculated by means of Eq. ( [39] l into Eq. ( |40l ). Then the given equation contains only the harmonics V'(O) and i/'(±2). Similarly, 
the transverse fields can be partially excluded from Eq. (|39l evaluating the harmonics ip{l) and ^(—1) by means of Eq. ( |40] |. In 
this case Eq. ( [39] l cannot be decoupled completely since it involves not only the longitudinal electric fields but also the higher 
harmonics of the transverse fields. Also it should be noted that the nonlinear response of the system is accompanied by the 
magnetic field generation according to Eq. ( flSt . It follows from this equation that the magnetic field is directed perpendicular to 
the plane containing the vectors k and Eg. 

In principle the dispersion equations of the perturbations can be deduced from Eqs. (|39] l and (l40l i by solving these equations 
by iteration to any order of accuracy. Taking, however, into account the smallness of the parameter ve/c, it suffices to restrict 
the analysis of the system ( [39b and ( |40l ) to the harmonics -E||(0), i?||(±2) and V'(O), '0(i2). In this case the longitudinal and 
transversal modes are decoupled and we obtain the following dispersion equations 



^\ik,-^) + i 
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'D±{k,uj — ujo) 'D±{k,uj + uJo) 
fc|| (fc, oj + vujo)ebii (fc, uj + I'ujo) 



fcf(fc,w)=0, (41) 



4c^ ^^ V\\[k,uj + jyuo) 

for longitudinal and transversal modes, respectively. Here ei,\\ {k, a;) = 1 + 6eb\\ (k, uj). 

A. Longitudinal modes 

Let us now investigate the dispersion equations (l4Tl i and (l42b in detail within fluid model (or cold plasma approximation) 
when the partial dielectric functions are given by 

SsAk,Lu) = &||(fc,a;) = 6s(uj) = -^. (43) 

Here tUp = Aimee^ /ra is the plasma frequency. Thus we consider only the high-frequency modes assuming that |a;| 3> fcwth, 
where wth is the thermal velocity of the electrons. We look for the solutions of the dispersion equations in the form lu ~ ijjr + ij, 
where Ur is the (real) frequency and 7 is the damping rate (when 7 < 0) or the growth rate (when 7 > 0) of the modes, 
respectively. In the absence of the laser field (ve — 0) the transverse modes are stable and their frequency is determined by 

^2 
ujl{k)=k^c^ + ^+iul (44) 

7& 

It is seen that the frequency of the ordinary transverse modes is modified by the charged particles beam effectively increasing the 
total plasma frequency of the beam-plasma system. Also it should be noted that the contribution of the beam to the dispersion 
relation of the transverse modes is ^ ^l/lb while for the longitudinal modes it is given by ^ '-^l/ll (see, e.g., Eq. ([3T|i). This 
is a consequence of the anisotropy of the effective electron mass with respect to the driving force acting either in longitudinal or 
transversal directions. The stability of the mode (04) can be easily understood taking into account the fact that the electric field 
vector in this mode is perpendicular to the beam and hence the work performed by this field on the beam particles is zero. 

The longitudinal two-stream modes are unstable in a long wavelength regime (see, e.g., Ref. {2^) with ^ fc ^ fee = u)c/uh, 
where lUc ~ Wp(l + S}^'^ /^bY^"^ and ^ = ujl/ui'i. Note that in practice ^ ^ 1 and Wc — '^p- Assuming that k <t^ kc the growth 
rate and the frequency of the two-stream modes read (cf. with Eqs. ( IA4b and ( IA5b ) 






7^S(fc) ^ fcu,^^^i^, a;;fS(fc)^_l^, (45) 



respectively. It is seen that the real frequency of the TS modes ~ kuf,, i.e., it is a frequency-locked oscillation, the frequency 
depending only the Cherenkov-type term and not on the natural frequency of the oscillations (^ ujp). The maximal value of 
7"""^ (k) is achieved at fc < fee 1281 and is given by 



/max „ , V '-' S 



24/3 



76 



(46) 



In the presence of the laser field (ve y^ 0) the high-frequency transversal and the low-frequency longitudinal modes, Eqs. (l44l l 
and ( l45T l. (l46T l, respectively, of the beam-plasma system are parametrically coupled according to Eqs. (HTt and (l42t . These 
equations can be satisfied only when one of the ordinary dispersion functions P|| (fc, w) or I?x(fc, w), becomes nearly equal to 
zero. This is not, however, sufficient to cause parametric TS instability, which occurs when at least two of the zeros of the 
dispersion functions merge, as in the case of the standard TS instability ll28ll38ll . In the case of the longitudinal waves there 
are three such situations: (i) Cherenkov-type coupling when w^ — kut,; (ii) D^^{k,uj) ~ and 'D±{k,uj ± uq) ~ 0; (iii) 
'D±{k, uj + cjo) — and I?±(fc, uj — ojo) — 0. The cases (i) and (ii) correspond to the resonant coupling, and the case (iii) to the 
nonresonant coupling. 

Consider the situation (i) which corresponds to the low-frequency (w ~ ci;Jf^(fc)), long-wavelength excitations. Close to the 

Cherenkov resonance, uj ~ fcu;,, the most important term in Eq. (ITU is involved in (5eb|| (k, w) (Eq. dsTl l). Consequently we look 
for the solution of the dispersion equation (l4ll in the form a; — ku\y + bj\, where |a;i| <C fcufe. In this case the instability occurs 
at kuh < Ws ~ Wp with the growth rate (wi — i-f) (cf. with Eq. ( l45T l) 

where ujs is the zero of the function F{lu), F{ujs) ~ 0, and 






1 1 



'D±{k,LL! - Uq) V±{k,LL! +UJ0) 



(48) 



It is noteworthy that the root fco of the transversal dispersion function, V^ (fc, oj — wq) — with uj = ku^, at wq < 2ujp may lie 
in the domain kub < uJs and the growth rate WT\ changes the slope close to this root. At fc > /cq the function F{ijj) becomes 
negative and the relation (l47T i is violated. However, at wq ^ Swp the formula (l47T i remains valid in the whole domain of the 
instability. Also the maximal growth rate is achieved at kub ~ LOg. Equation ( |47] | is clearly invalid in this case and more rigorous 
treatment of the dispersion equation yields 



7max — „4^3 



V3 r i 20.2 1 1/3 



7,' |f 'K)l 



(49) 



Here the prime indicates the derivative of the function with respect to the argument. 

We have considered above the low-frequency and long wavelength regime when the dispersion properties of the system are 
strongly determined by the beam characteristics (density and the energy). Next we consider high-frequency {uj > ujJ^ (fc)), short 
wavelength regime with kui, > uig. This case corresponds to the situations (ii) and (iii) introduced above. Assuming low density 
beam with rib ^ n^, we note that the role of the beam in the dispersion properties of the system becomes less pronounced in this 
high-frequency regime, and as a first approximation the terms proportional to cj^ can be neglected in Eq. dTIT l. This is a regime 
of purely parametric excitations in a plasma. 

The resonant coupling (ii) occurs when the following resonance condition is satisfied: 

Wo — Wp + uj^{k). (50) 

In either case, \uj\ is assumed to be much smaller than wo and cj^(fc). We can them make the resonance approximation for the 
dispersion function I?_l(A:, oj ± wq) to obtain 

2?i(fc,c^±c^o)-T;^(2wo-^)(^±'^), (51) 

where 5{k) = ujq—uj± (k) is the mismatch of the laser frequency from the frequency of the natural transversal mode. Substituting 
Eq. ( BTl i into (|4TI) and neglecting the contribution of the beam in the dispersion relation of the left hand side of Eq. (l4ll yields a 
cubic equation for w^, 

-'(-'-"?)("' -*^) + ^;^-° (52, 
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which can be solved analytically. However, we restrict ourself to the simple and qualitative solutions of the dispersion equation 
and more rigorous numerical calculations will be presented in Sec.lVll 

As mentioned above there are two types of solutions of the high-frequency (uj > W|F^(fc)) dispersion equation (|52] |. First 
we consider the resonant solution assuming that 5 ~ Wp > 0. Note that this type of instability resembles the resonant decay 
instability. The sum frequency of the excited modes is exactly equal to the laser frequency uo (see Eq. (fSOll). Introducing the 
frequency mismatch A(fc) — S{k) — Wp (A <C ujp, 5) the dispersion equation ( |52] i for the resonant growth rate yields 



(fc) = iVe^fc^c2-A2(fc), (53) 



7,.. 2 



where e = {ve/'2,c){t — 1/2) ^/^^ -j- = a;o/ajp. The resonant unstable mode exists only at r > 2. The maximal growth rate 
7r,max IS achieved sik — fc^.max with 



^^ '^ '-, (54) 



C 1 — e 



Wp / r(T-2) + e2 

7r,max = yey ^3^5 • (55) 

Now in this high-frequency domain the characteristics of the instability are only weakly sensitive to the beam density and 
energy (7;,) but essentially depend on the laser intensity and the frequency. It is seen that the maximal growth rate is scaled as 
(at T ^ 1) 7.r,max ^ [-^L ('^p/'^o)] ^^^^ whcrc II is the RF intensity. It is also noteworthy that at r ^ 1 the position fc^.max of 
the maximum of the resonant growth rate is independent on the laser intensity (fcr,max — ^o/c). Equations (l53Tl-(l55Tl can be 
compared with the growth rate of the ordinary two-stream instability, Eq. ( |46] l. Assuming, for simplicity, high-frequency laser 
field ((jJq ^ (jjp) one obtains that at ve/c > \/S/T{4:S,/^^y^'^ the growth rate 7r,max exceeds j^^- Note that the last inequality 
is easily fulfilled for the REB. 

In the nonresonant case (iii) assuming that |a;| ^ cUp from Eq. (|52] | one obtains a quadratic equation for w^. A simple 
analysis of this equation shows that in this case the instability occurs at ~Sm ^ S ^ ujo — ujp, where Sm = ujp[2ei{kc/ujp)]^/'^ 
with ei — {ve/'2c){t — 6/2ujp)^^^^. Two distinct branches of the instability should be considered separately. At the positive 
frequency mismatch, ^ 5 ^ wo — ^p, (or k ^ k2 = {u}p/c){T^ — 1)^/^) the instability is almost aperiodic (i.e. iVr — 0) with 
the growth rate 

Al/4 , , , s 1/2 

^='—i^,/pTSl-0') , (56) 

while at the negative values with —Sm ^ S ^ (k2 ^ fc ^ fci, where 5 = —5„i at k = ki) the instability is periodic (i.e. 
iOr 7^ 0) and the growth rate becomes 

., l^r^V.3/2 ,.,3Ml/2 



{sT - 1^1^/') . (57) 



The real frequency of the unstable mode (|57T l is obtained by changing the minus sign in this formula by plus sign. We note that 
the latter nonresonant case with ^ < resembles the oscillating two-stream instabiUty. 

Thus summarizing this section we emphasize that in the spectrum of the longitudinal unstable modes there are basically three 
domains with strictly different properties. The "long wavelength" domain with k < cjp/ub (we denote this parameter regime 
as Domain I) is basically determined by the beam density and the energy (7f,) and corresponds to the TS instability (see, e.g., 
Eqs. (|46] l and (|47|i). The intermediate (Domain 11, k ^ ^2) and the short wavelength (Domain 111, fc2 ^ fc ^ fci) regimes 
mainly depend on the laser intensity and are only weakly sensitive to the beam parameters. In 11 and 111 the growth rates can 
be approximated by Eqs. ( |56] | and (l57T i. respectively. In addition in Domain 11 at luq > 2ujp it is possible to witness a resonant 
excitation of the unstable modes with maximal growth rate (|55] l which resembles the resonant decay instability. Finally, the 
Domain I may merge to Domain 11 at fc ~ ojp/ui, while the Domain 11 merges to Domain 111 at zero frequency mismatch, 6 = 

(k = fc2). 

B. Transversal modes 

In this section we turn to the investigation of the unstable transversal modes generated in the beam-plasma system. Our 
starting point is the dispersion equation (l42l i for these modes. As in Sec. lIV Al we adopt here a cold-fluid approximation when 
the partial dielectric functions of the beam and the plasma are given by Eqs. (ISTT l and ( |43] l, respectively. An inspection of 
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the dispersion equation (l42T i shows that there is only one resonant coupling between different modes. This is the situation when 
'D±{k,uj) ~ andl?||(fc, w — wg) ~ 0. This system of equations require that w ~ uJo+uJs{k) + Au!{k) andijj±{k) ~ uJo+ujs{k), 
where Aw <C wq + ujs, i^± (k) is the frequency of the ordinary transversal modes, Eq. (l44l . and tOs is the real root of the ordinary 
dispersion equation, 'D^^{k,LUs) — 0, for the longitudinal modes. Note that the resonant condition uj± ~ uo + ojs cannot be 
satisfied in the domain where the two-stream instability occurs, where Wg (/c) = wjf ^ (fc) + ^7^^ {k) (see, e.g., Eqs. (|45] | and (|46] |) 
is a complex quantity. Therefore it is expected that the resonance occurs at short wavelengths (fc > kc) and at high-frequencies 
(u! > a;jf^(fc)). Inserting a; ~ ajp + Wg + Auj and uj± ~ wq + i^s into Eq. ( |42] | and neglecting a small term depending on the 
frequency uj + ujoin the right-hand side of this equation, for the maximal growth rate we obtain 




In addition, the resonant instability occurs only for a negative derivative, ■^'D\\{k,ujs) < 0, of the longitudinal dispersion 
function. An analysis shows that there is only one real root Ug = lj^-^ of the equation I?||(fc,Ws) = which satisfies this 
condition (see Appendix lAl for details). This root is negative and is represented here as uj~i{k) = -~ujpg{k), where the function 
g{k) is positive and decreases monotonically from (1+^/7^)^^^ (at fc ~ 0)to 1 (atfc ^ kc), see Eqs. (IA6I) and dAlb . respectivelv. 
The maximal growth rate (ISST i is reached at A; = fcmax which is determined by the resonant condition above, uj±{k) = 
Wq — ^pg{k)- Assuming short wavelengths (fc > k^) this equation can be solved iteratively. The result reads as 

k'^ r^ ^ f 

^H^^[r-5M]^-^-l, (59) 

where k^c^ = a;p[r(T — 2) — C/7b]- Note that the resonant instability occurs only at high-frequencies of the RF, t > Tc = 
g{K) + (1 + ihhY^'^ (for a low-density electron beam this condition is roughly equivalent to ujq > 2ujp). 
Substituting Eq. i5% into Eq. ( |58] l we arrive at 



7r,max ^, Ve [t - g{K)f - ^/"fb "1 1 



ujp 4cV gin)[T-giK)] l + [^l^l)H{n) 

with 



JTTTT^ (60) 



Hin) 



KUb/iOp + g{K) 



(61) 



The maximal growth rate for the resonant instability is strongly simplified for a very low density {n\j <C rie) or for an 
ultrarelativistic electron beam (7^ ^ 1). In this case the instability occurs at r > 2 and /cmaxC — u)p[t{t — 2)]^/^. The growth 
rate 7r.max is then simply given by 



7r,max Ve T (t - 2) 

u)p 4c V r — 1 

and is completely independent on the beam parameters (this is a regime of a purely parametric instability). Note that similar 
to Eq. ( fSST i for the longitudinal resonant unstable mode the resonant growth rate (l62T i at wq ^ LOp is scaled as 7r,max ^ 
[/^(wp/wo)]^^^- Thus the growth rates of the resonant longitudinal and transversal unstable modes may be of the same order 

Besides the resonant mode there also exist two nonresonant transversal modes. The dispersion equations for these modes 
follow fromEq. ( |42] | and aregivenby 2?||(fc,a;±a;o) — which implies that a; = ^pajo+'^s + Aa;± (with \i-u±\ ^ | =Fajo+'^s|)- 
Here ujs = ujJ^ + 17^^ (see, e.g., Eqs. ( l45T l and ( IA4l i. (IA5l l) is the solution of the longitudinal dispersion equation in the domain 
of the two-stream instability which occurs at the low-frequencies (a; ~ ujT^{k)) and at the long wavelengths (fc ^ kc). 

The quantity Aa;± can be roughly estimated using the dispersion equation (l42T l. First, in the leading order we represent 
the dispersion function in the form 2?||( fc, lo ± loq) ~ Aaj±^2?|| (fc, lOs). Second, the derivative of the dispersion function is 
estimated employing Eqs. (IA4b and dASI l. Assuming a low-density (ri;, <C Uc) or an ultrarelativistic (7;, ^ 1) electron beam this 
quantity in the leading order of the dimensionless parameter ^/7i^ is represented as -§^T>\\{k,ijJs) ~ {2i/i^^){LOp/LoJ^Y. Then 

„3 , 



inserting these results into dispersion equation (142) in the leading order of 1^/7^ one finally arrives at 

8 (a;f)2[(c.oTwf)2-c.i 



AU± ~ -I ^ , ,TRN9r/ TSN9 -7T- (63) 
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The real frequencies and the growth rates of the nonresonant modes are simply given by ujr — T^o + wjf^ and 7 = 7^^ + 

Im[ALL'±], respectively. It should be emphasized that the frequency shift ( l63T l is valid when Im[Aa;±] <C 7^^. Moreover, since 
the nonresonant unstable modes appear in the domain of the two-stream instability with fc ^ fc^ the frequencies JJ^ and uj±_ 
in the brackets of the denominator of Eq. (l63T l can be neglected and hence Aaj_|_ ~ Aw_. The expression (|63] | can be further 
simplified recalling that luJ^ ~ kub and the ratio 7^^/ajjf^ ~ {^HbY^"^ is almost independent on k (see, e.g., Eqs. (IA4l l and 
(IA5b ). Therefore the growth rate of the nonresonant modes is given by 



7(fc)^7^^(fc) 



l_i f — V f^ 



8 \Ub J \ujo 



(64) 



It is noteworthy that this growth rate is proportional to the growth rate 7^^ (fc) of the standard two-stream instability and only the 
factor in the brackets depends on the laser intensity. Thus 7(fc) is weakly sensitive to the laser intensity (because in the present 
approximation ve <C Ub and uiq > uip) and is mainly determined by the beam-plasma interaction. In addition, the nonresonant 
unstable modes do not disappear with decreasing the laser intensity as it occurs for the resonant ones. 

V. LONGITUDINAL POLARIZATION OF THE RF (Eo | k) 

With the theoretical formalism presented in Sees. |II]and|III] we now take up another configuration of the laser polarization. 
In the following we study in detail the parametric two-stream instabilities in the laser irradiated plasma when the polarization 
vector of the laser field is parallel to the wave vector k of the excitations (Eq || k) assuming again that the beam is directed in 
the direction of k. It is expected that the beam-plasma and the laser-plasma unstable modes are strongly coupled in this regime 
compared to the transversal configuration (Eg J- k) since the electrons are effectively driven by a laser radiation in the direction 
of k (II Uf,) in this case. 

Our starting point is the general equation (l25t for the harmonics which for the configuration Eq || k || Uf, is decoupled into 
two independent equations for the longitudinal (_E||) and transversal (Ej^) electric fields, respectively, 

CO 

^||"^(0)+ E £;|f^(0)J,_,(C)J„-.(C)&||(s) = -<J„o^l|(0)fct||(0), (65) 



£,s— — oo 

A;"^c^ — cj^ 7b 

Here ( = a(k • e) = ±ka and the other notations have been introduced in Sees. |II] and |III1 As in the preceding Sections we 
consider throughout an electron plasma neglecting the dynamics of plasma ions. 

To exclude the electric field harmonics En^''' in Eq. ( |65] ) we multiply both sides of this equation by Jn-p{C)Je-piO ^nd 
perform a summation over n and p. As a result we arrive at 

i?j")(o) = -£;|i(o)fc,||(o) £ t+feTir ^^^^ 

Thus the harmonics ^Ey are completely expressed by the total field E\\. Also in deriving Eq. < |67] > we have used the summation 
formula for the Bessel functions [42] 

Y, Js^e{C)Js-n{0 ^ Sin- (68) 

s— — oo 

Next, in Eq. ( l67] t the frequency uj is replaced by cu — nujo and using Eq. (l24] t we perform a summation over n. This yields an 
equation for the amplitude £"11 (0) of the electric field 

CSO 

[l + fcb||(0)]i?||(0)= Yl ^l|W*KC)fcfc||M), (69) 
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where 



^^(0= E 



fc||(-s) 

l + Seiii-s) 



Js-iiOJsiO- 



(70) 



It is noteworthy that in the Eo || k geometry the system of equations (l65T l and (l66T l involves all (infinite number) the harmonics 
of the electric field whereas in the transversal case, Eq -L k, considered in Sec.|lV]each harmonic E^'^'' (0) with given frequency 
uj connects only to the nearest neighbors _E("^^'(0) and E'-''^^^\0). These features are the peculiarities of the specific laser 
polarization directed parallel or perpendicular to the wave vector k. 

Similarly, it is possible to exclude the harmonics E]" of the transversal electric field from Eq. (l6&i and derive an equation for 
the electric field E^, 



Ei(0) [k'c' + ^-u'] =uo 



b_ 

lb 



E 



E 



-(■) 



{u + sa;o)^ 



Jt-s{0J-s{05e^{s) 



(71) 



First, we note that at vanishing laser intensity from Eqs. (I69]l-(l7ni we recover the standard dispersion equations for the longitu- 
dinal, Pji (/c, oj) = 0, and transversal, I?^ (fc, lo) = 0, waves given by Eqs. ( l32b and ( |38] ). respectively. Second, within cold-fluid 
approximation (the dielectric functions are determined by Eq. ( |43] l) using the summation formula ( |68] l for the transversal modes 
we arrive at the same dispersion relation as in Eq. (04). Therefore within this approximation the laser field has no influence on 
the dispersion properties of the transversal modes which are stable in this case. Third, in contrast to the longitudinal modes ( |69l ) 
the transversal ones are only weakly sensitive to the beam parameters (see left-hand side of Eq. (jTT)) ^nd the instability of these 
modes is parametric in nature. These instabilities have been studied previously in Ref. |44]. Consequently, in the following we 
consider throughout only the dynamics of the longitudinal modes E\\ . 

As was mentioned in the previous sections the dispersion equation of the perturbations can be deduced, in principle, from 
Eq. ( |69] l by solving this system of equations by iteration to any order of accuracy. However, to gain more insight we consider 
the long wavelength limit of Eq. ( |69] l when the parameter \C,\ = kais small, which suffices to restrict the analysis of the system 
I to the harmonics E^\ (0) and E\\ (±1). In this case we obtain the following dispersion equations for the longitudinal modes 



Vii{k,Lo) 



k'^a' 



5eb\\{k,uj)[Pi{k,Lo) + P-i{k,uj)] . 



(72) 



where 



F±i(fc,w) 



(5£b||(fc,i:^±i:jo) +£||(fc,cj) 
'D\\{k,ijj ± Wo) 



e||(fc,w±a;o) 



e\\{k,uj) 



(73) 



The dispersion equation ( |72] i can be compared with Eqs. (ITTT l and (l42b obtained for the Eq J- k geometry. It should be 
emphasized that unlike the Eq J- k geometry the dispersion equation (|72] | (see also the more correct relation ( |69) ) in the absence 
of the charged particles beam (fc^y — 0) yields the standard dispersion equation e|| (fc, lu) — for the longitudinal modes. Thus 
in this case the laser radiation with Eq || k does not influence the dispersion properties of the plasma in the k direction but it 
affects the dispersion relation in the transversal direction (see Eq. (ItTI)). This result sounds paradoxical, considering that the 
plasma oscillations should be effectively driven by the laser field in the Eg || k configuration. However, let us recall that the 
dynamics of the plasma ions is completely neglected here. In reality, the laser radiation with Eq || k stimulates low-frequency 
(typically with the ion plasma frequency) electron-ion coupled oscillations 1.37^ .38.1 and the dispersion relation is not simply 
given by the equation ey (fc, w) = 0. 

The simplest way to investigate the parametric TS instabilities determined by Eq. (l72T l is the cold-fluid model when the 
dielectric functions are given by Eq. (|43) . In this specific case the function '^i{C) is evaluated analytically in Appendix IB] 
(Eqs. ( |B3l l and ( IB4I )) using the Newberger's summation formula ll46ll . Since the function ^^ (C) decays exponentially with (. the 
harmonics E\\ {€) in Eq. ( |69] l and the corresponding dispersion relations can be effectively evaluated numerically to any order of 
I. Furthermore, within cold-fluid approximation it is possible to derive a dynamical equation for the complex amplitude of the 
excited waves. This is done in Appendix ICl see Eq. ( IC6I ). 

One characteristic feature of the dispersion equation (|72i for the longitudinal modes in Eq | | k geometry is the absence of the 
contribution of the transverse modes. We consider first the situation (i) (see Sec. lIV Al l which corresponds to the low-frequency 

Close to the Cherenkov resonance, lo ~ fcufc, we look for the solution of the 



(a; ~ a;jF^(fc)), long-wavelength excitations 
dispersion equation (l72T l in the form cj = ku\j + cji, 
growth rate (uji ~ ij) (cf. with Eqs. ( |45] l and ( |47l )) 



where \uji\ ^ kui,. In this case the instability occurs at kui, 
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ujp with the 
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where 
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[Pi(fc,c^)+P_i(fc,c^)] 
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(75) 



Note that Eq. ( l74l ) is vahd far from the roots of the equation e{uj ± cjq) = '^b/lb^O' where e{Lu) = 1 + 5e{uj) is the dielectric 
function in a cold-fluid approximation. 

The maximal growth rate is achieved at kuf, ~ ujp. Similar to Eq. (l47T i. the relation (|74] | is clearly invalid in this case and more 
rigorous treatment of the dispersion equation for the frequency correction ui yields a fourth order algebraic equation 



Here 
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UJl 



UJp$. , 

:2 7^ZiiPi 



8ui 7,^r4 



e(wo ± Ldp) 



P^i) 



:{ujo±ujp)- ujlhl^l ' 



(76) 



(77) 



and the laser frequency ujq is not too close to the value 2ujp[l + (^/27J^)]. It can be shown that the second term in the right-hand 
side of Eq. ( |76] ) is systematically smaller than the first one. Neglecting this term we arrive at the maximal growth rate 'j'^^^ of 
the two-stream instability, see Eq. ( |46] |. Thus as expected in the situation (i) the maximal growth rate is only weakly affected by 
theRF. 

Next we consider high-frequency (a; > ujj^ (fc)), short wavelength regime with kub > Wc- The nonresonant coupling similar 
to (iii) is determined by the intersection of the different roots of the dispersion equations 2?|| {k, to ± ujq) ~ 0. If ujsi and uJs2 
are two different roots of the parallel dispersion function, I?|| (fc, uJsi) = I^n {k, lUs2) = 0, the above mentioned intersection 
of the roots simply yields cugi ~ 0Js2 + 2wo. This equation implies that both roots ujsi and uJs2 should be real which is only 
possible outside the domain of the two-stream instability {kub > ujc^ see also AppendixlAli. In addition, the frequency shift Acj 
determined by the relation uj + ujq = LUgi + Alu (or alternatively cu — luq = uj^^ + Aw), is also real and can be calculated from 
the dispersion equation (|72] |. Therefore in this case the nonresonant coupling merely shifts the real frequencies of the modes and 
does not cause any instability in the beam-plasma system. 

Instability may occur in situation (ii) with the resonant coupling. In this case Py (A;, w) ~ and 2?|[ (/c, uj ± wo) — 0. Again 
introducing two different real roots los\ and los2 of the parallel dispersion function we consider the frequency shift Aw with 
u) ~ ujs2 + Aw, such that |Aw| ^ |wsi,2|- The resonant coupling occurs when uJs2 = <^si ± i^o- Substituting this relation and 
the frequency ut into dispersion equation (|72] | and taking into account the smallness of the quantity Aw we obtain 



Aw^ 



fc2a2 



[s{uJs2) ~ £{^sl)Y 



£l?||(fc,W,i)£l?||(fc,W,2) 



(78) 



It is seen that this expression is symmetric with respect to the exchange of the roots Wsi and ujs2 and yields an unstable mode 
if the derivatives of the longitudinal dispersion functions in the denominator of Eq. JTSl l have different signs. The resonant 
coupling condition, ujs2 = ^si ± ^o. and the restriction on the signs of the derivatives of the longitudinal dispersion functions 
in Eq. dTSl l reduce the possible candidates for the quantities w^i and 0Js2- From Appendix lAl it follows that there are only three 
choices; a) Wsi(fc) = w^^(fc) and uJs2{k) — oJrzik) with Wr3(A:) = wq + w^^(A:);b) Wsi(fc) — w^^(fc) and Ws2(A:) = (jJr2{k) with 
^r2{k) = Wo + uj^i{k); c) Wsi(fc) = '^.^ilfc) and ujs2{k) — i^ri(^) with co^i{k) — ujq + w^j^(A:). In these cases the derivatives in 
the denominator of Eq. ( iTST i have different signs and the corresponding modes are unstable. The positions fc^ax of the maximal 
growth rates of these modes can be determined from the resonant coupling condition. Using the asymptotic behavior of the roots 
ujfi (k), ujr2 (k), and w,.3 (k) at fc ^ kc, Eqs. ( IA1I )-( |A3] ), respectively, we introduce four new functions f± (fc) and h± (k), which 
are determined through relations 



0Jr2{k) 
UJrz{k) 



kubT ^p^n^h±{k), 



ujfiik) = ±Wp 



7& 



"P 3/2' 

7fc 
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(79) 
(80) 



Then in the cases a) and b) with Ws2 = Wgi + wq the resonant coupling conditions for the determination of /cmax in a leading 
order of the parameter £,/j^ yields a pair of the transcendental equations 



kub = Wo ± Wp ± —^u}ph±{k) 
lb 



(81) 
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which can be solved iteratively. Here the minus and plus signs are related to the cases a) and b), respectively. Since the functions 
h±{k) at fc ^ fcc behave as h±{k) = 1 + 0{{ujp/kui,)'^) (see Appendix lAl) within zero order the last term in Eq. (IHTj can be 
neglected which yields k± = (wq ± ujp)/ub- Substituting this value into the arguments of the functions h±{k) in the last term 
of Eq. dSTT l one obtains the corrections to k±. The maximal growth rates are obtained from Eq. (fTFt . where k = k^-aax- In the 
leading order of £,/^l the result reads 



-V^ „ / t \ 1/4 

'r,max ^E I S 



4W6 \% 



T±2 



r±l 



[h±{n±)f\ (82) 



For an estimate the approximate expressions /i±(k±) ~ 1 + (1/2) (r ± 1) ^ for the functions h±(k) can be used. Here 7^ ^^^ 
and 7,tjjja^x are related to the maximal growth rates in the regimes a) and b), respectively. Let is note that for a validity of 
Eq. (|82] | in the regime a) the laser frequency cjq should not be too close to the plasma frequency. From Eqs. (|72]) and ( f78b it is 
straightforward to obtain the profiles of the resonant growth rates in the regimes a) and b). Introducing the frequency mismatch 
6 = ujs2 — Wsi — Wo these profiles are determined by 



itik) - y4(7±nax)'-<52(fc). 



(83) 



Here 7^max ^'"^ ^^^ maximal growth rates (see Eq. (l82t ) which are achieved at (5 = 0. It is seen that the quantities 7^ (fc) and 
7+(fc) vanish at ki-2 — k- ± 2j^^^^^^/ub and k3-4 ~ k+ ± 2^^^^-^^Jub, respectively. 

Consider now the regime c) when ^^2 = oJsi — wq. In this case the resonant coupling condition reads 

Wo-2wp=Wp43[/+W + /-W]- (84) 

It is clear that this relation can be satisfied only at wq > 2ujp since the functions f±{k) are positive (see Appendix lAJi. On the 
other hand wq should be sufficiently close to 2Lu,p because the dimensionless parameter £,/jb i^ small. Assuming for instance 
T— 2 ^ '?/7b and employing the results of the Appendix|A]the solution of Eq. (|84]l reads as fco^b/'^p — (S./lbV^^i'''^'^)^^^^ ^ 
1. This value of k determines the position of the maximal growth rate in c) which is obtained from Eq. (iTST i and is given by 

7r,max Ve T-2 ^^ Ve ^ 

Up Ub T Ub T% 

Confronting this relation with the growth rates 7i^,jiax we conclude that 7r,max <C 7^max- I" addition it should be noted that in 
general fcg 3> n± and the growth rate 7r,max may be strongly shifted towards very large k values. 

Unlike Eq J- k geometry the beam-plasma and the laser-plasma parametric modes are strongly coupled here. As a result the 
resonant modes (|82] | and dSSl l depend essentially on the laser intensity and the beam density. 

Let us now compare the growth rates in the regimes a) and b) for the resonant unstable modes with the quantity 7,^ax assuming, 
for simplicity t ^ 1 (ojo ^ Wp). In this case 7^,„ax — Ti^max- It is seen that 7^niax exceed the growth rate 7jjax of the two- 
stream instability at sufficiently intense RF, ve/c > 3.02(7^ax/'^p)^^^- It is clear that this condition requires very low density 
beams (compared to n^) and is increasingly difficult to fulfill with increasing Ub- Also the growth rates in Eq. (|82] | should be 
compared with the growth rates of the resonant longitudinal (Eq. dSSl l) and the transversal (Eqs. (|60] | and (l62t ) modes as well 
as with the nonresonant longitudinal modes (Eqs. (|56] | and ( fSTJ l) excited in Eq J- k geometry. Again assuming, for simplicity 
the Wo ^ ujp limit we conclude that the unstable modes grow much faster in this geometry, where the parametric effects are 
more pronounced. However, it should be emphasized that the resonant unstable modes in Eq J- k geometry are only effectively 
excited starting with the threshold frequency ujq ~ 2ujp of the laser radiation. Below this threshold (with luq < 2wp) only the 
unstable mode ( |64l ) in the long wavelength domain and the modes ( |56] ), ( |57] l, and ( [82] l in the short wavelength domain are excited 
in Eq J- k and Eq || k configurations, respectively. 

VI. NUMERICAL TREATMENT 

Using the theoretical findings of Sees. UllWl we present here the results of our numerical calculations of the growth rates for 
the longitudinal and the transversal unstable modes assuming the transverse (Eg J- k) and the parallel (Eq || k) configuration 
of the RF amplitude Eq with respect to the wave vector k. The calculations have been done for an electron beam with a small 
dimensionless density parameters £, = (ujb/'^p)'^ = rib/ue — 0.1 and <^ — 0.3 and for a relativistic factor jb = 5. For the 
laser intensity parameter a = v^/c^ we have adopted the values a = 0, 0.01, 0.1, and 0.2. It is convenient to represent the 
laser intensity parameter in the form a = {Il/Iq)Xq, where Iq = 1.37 x 10^^ W/im^/cm^ and the wavelength (Ao) and the 
intensity (7^) of the laser field are measured in units /im and W/cm-^, respectively. The laser frequency is measured in the units 
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7^ = 5; n|^ = 0.1n^; u^= 1.2(t 
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FIG. 1: (Color online) The growth rate 7|| (in units of plasma frequency cop) of the longitudinal unstable modes in terms of kub/i^Jp obtained 
by numerical solution of the dispersion equation (I41t in Eo -L k configuration for 76 = 5, nt = O.lrie, ojq ~ l-2ujp, a = (solid line), 
a = 0.01 (dashed line), a = 0.1 (dotted line), a — 0.2 (dash-dotted line). The left and the right panels correspond to the branches 1, II and 
m, IV, respectively, introduced in the text. 



of the plasma frequency, t = uJo/utp > 1. In our numerical calculations this parameter varies in a wide interval, 1.2 ^ t ^ 4. 
Throughout in this section the growth rates are measured in units of plasma frequency cUp and are calculated as a function of 
kub/ojp for several laser intensities and frequencies. Note that the chosen parameters both for electron beam and RF are typical 
for FIS for inertial confinement fusion [7]. Assuming, for instance, radiation field with Aq — 0.5 /im, the parameter a ~ 0.2 
corresponds to the intensity II — 10^^ W/cm^ of the RF. 

First we consider the transverse geometry with Eo J- k. In this case the basic properties of the dispersion relations for the 
beam-plasma system have been studied in Sec.|IV] In general we have found that the simplified treatments of Sees. |IV| and FVl 
agree qualitatively well with the exact numerical solutions. However, it is clear that these simplified treatments are not capable 
to resolve all details and the branches of the spectrum of the unstable modes in the cj-fc plane. 

Within Eo J- k geometry we now consider the case of the longitudinal unstable modes (7 — 7|| ) when the dispersion relations 
are determined by Eq. (l4ll . As mentioned in Sec. II V A I this is a regime when the purely two-stream and the parametric modes 
are only weakly coupled and their spectra are well separated in the uj-k plane. Therefore the different modes (beam-plasma 
or parametric) are only weakly sensitive either to the electron beam parameters or the intensity of the RF. From Eq. (HTi it 
is seen (see also the simplified version of this relation, Eq. ( |52l )) that in a cold-fluid approximation there are ten solutions of 
this equation, but only some of them correspond to the unstable modes with 7 > 0. In addition, these unstable modes are 
excited with different real frequencies. To demonstrate this feature in Fig.[T]the growth rates are shown for the laser frequency 
ujQ — l.2ujp. In this case we have found numerically that there are only three solutions which coiTespond to the unstable modes, 
and as an example two solutions are shown in the left and right panels of Fig. [1] The different curves correspond to the laser 
dimensionless intensities a ~ 0.01 (dashed lines), a = 0.1 (dotted lines), and a = 0.2 (dash-dotted lines). The solid lines with 
a — represent the growth rate of the standard two-stream instability, Eq. ( |45] ). Two panels of Fig.[T]correspond to the modes 
with different real frequencies w,.. We denote tentatively the solutions with different ujr as the branches I, II, III etc. In the left 
panel of Fig.[T]up to the value k2c/u]p ~ (r^ — 1)^/^ ~ 0.6 {k — fc2 corresponds to the vanishing frequency mismatch, S — 0, 
introduced in Sec. |IV At the real frequency is a;^ = (branch I), while at A;2 ^ fc ^ fci it is given by cj^ — ojg{k) (branch II). 
In Fig.[T]the growth rates sharply tend to zero at k = ki. [For the approximate definition of the quantity fci see paragraph above 
Eq. (fSTbl. The spectrum ujg{k) corresponds to the real frequency of the nonresonant longitudinal modes deiived approximately 
in Sec. lIV Al for the negative frequency mismatch (6 < 0). In the approximate form it is given by Eq. (|57i , where, however, the 
minus sign has to be replaced by the plus sign. In the approximate treatment of Sec. IIVAI Eqs. (|56] | and (ISTT i correspond to the 
branches I and II, respectively. Thus the boundary between the branches I and II is determined hy S — (or k — k2). In the 
branch I the growth rate at fc < A:2 increases almost linearly with fc, j{k) ~ A;i;£;[2(r^ — 1)]^^/^ in agreement with Eq. ( |56] l. 
Finally, in Fig.[T](right panel) the real frequencies of the modes at ^ /c ^ fee and kc ^ k ^ max[A;c; ki] coincide with the real 
frequency of the standard two-stream instability, utr — u!j^{k) (branch III, see Eq. ( |45] |) and ojr = ojq + (^±{k) (branch IV), 
respectively. Here uj±{k) is given by Eq. ( |44| ). Let us note that the boundary k = k2 between I and II does not depend on 
the intensity of the RF while the upper boundary k = ki of the branches II and IV (at fci > kc) is shifted towards shorter 
wavelengths roughly as ki ~ fc2[l + 0(q;^/"^)] with increasing laser intensity. It is clear that at smaller intensity of the RF when 
ki < kc the branch IV disappears. As a general rule we observe that at kc ^ k ^ ki the growth rates in the branch IV are the 
part of the branch II, except the case of the low intensity RF with ki < kc (cf. the dashed curves in Fig.[T]with a — 0.01). In this 
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FIG. 2: (Color online) Same as in Fig. [T]but for the higher laser frequencies ujq = 2ujp (left panel) and cjq = 4tjp (right panel). Note the 
different scales in the left and the right panels. 
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FIG. 3: (Color online) The growth rate 7x (in units of plasma frequency ujp) of the transversal unstable modes in terms of kut/uip obtained 
by numerical solution of the dispersion equation ( I42t in Eo ^ k configuration for 7;, = 5, rit = O.lrie, loq = 2up (left panel), uio = Auip 
(right panel), and f or q = (solid line), a — 0.01 (dashed line), a = 0.1 (dotted line), a = 0.2 (dash-dotted line). Note the different scales 
in the left and the right panels. 



low-intensity limit the branch IV disappears while the branch III is nearly the same as the standard branch for the two-stream 
instability (solid curves in Fig. [U. Thus, at ujq > tOp and at small intensities of the RF (fci < fee) the parametric two-stream 
instability occurs in the branch III with the growth rate ~ 7^^(fc) which is only weakly affected by the RF. At higher intensities 
of the RF (with fci > kc) a new unstable branch IV is formed. The branches I, II and IV are formed due to the parametric 
excitations and are almost insensitive to the electron beam. 

Next in Fig. |2]the longitudinal growth rates are shown for the higher laser frequencies ujq = 2ujp (left panel) and wq = 4ajp 
(right panel). As expected the different branches shown in Fig.[T]are mixed with increasing wq and, in addition, more and more 
new branches for the unstable modes are formed. As an example in Fig.|2]only two solutions of the dispersion equation (ITTT i are 
shown which involve the basic features of the branches I, II, III, and IV introduced above. Note that in these particular examples 
with higher laser frequencies, ^2 — ij-iq/c exceeds the upper boundary kc of the two-stream instability, and therefore fci > kc for 
an arbitrary intensity of the RF. As has been pointed out above the approximate growth rates given by Eqs. ( |56t and (|57t are not 
capable in these regimes to resolve all specific branches shown in Fig.|2l The most important feature shown in Fig.|2]is that the 
domain of the instability in the fc-space is broadened accompanied by an increase of the maximal growth rate with increasing 
cjQ- It is also noteworthy the formation of the resonant growth rate in the spectrum of the unstable modes at ujq > 2ujp derived 
approximately in Sec. IIV Al see Eqs. (|53]|-(|55]|. Let us recall that the resonant unstable mode is not excited at cjq < 2ujp (see 
Sec. IIV Al l and hence this mode is not visible on Fig. [l] As predicted by Eqs. (|53])-(|55]) the resonant coupling at luq — 2ujp 
(Fig. 121 left panel) is only weakly pronounced with the maximal growth rate 7r,max/wp ~ a/ 6 and /cr.max — ('^p/c)(a/3)^/^ 
while at higher frequency ujq = Aojp (Fig.|2] right panel) it is strongly increased and is shifted towards the short wavelengths. 
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FIG. 4: (Color online) The growth rate 7|| (in units of plasma frequency cop) of the longitudinal unstable modes in terms of kut/(jJp obtained 
by numerical solution of the dispersion equation ( 172) with Eq. l l73t in Eo |[ k configuration for 7^ = 5,nb — O.lne, ujo — 1.2LOp (left panel), 
ojo = '2ijJp (right panel), and for a = (solid line), a = 0.01 (dashed line), a = 0.1 (dotted line), a = 0.2 (dash-dotted line). 
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FIG. 5: (Color online) Same as in Fig.|4]but for tJo = 4:iUp, nt = O.lrie (left panel) and nt — O.Srie (right panel). 



7r,max/i^p ~ (ar/4)^/^ with fcr.max ~ "^o/c- In Fig. |2](right panel) the resonant growth rate is determined by the curve with 
the maximum located around kub/ujp ^ 3. 

The growth rates for the transversal unstable modes are demonstrated in Fig.[3]for the Eg -L k geometry and for cjq = 2wp (left 
panel) and cjq = 4a;p (light panel). These growth rates are obtained by the numerical solution of the dispersion equation (|42] | 
for the transversal modes. The results for the smaller laser frequencies (with cjq < 2wp) are not shown in Fig. [3] This is 
because only the nonresonant modes with the growth rates approximately given by Eqs. (|63] | and (|64] | are possible in this case as 
discussed in Sec. jIVBI We have found numerically that in this frequency regime the growth rate only weakly deviates from the 
growth rate 7^^(fc) (see Eq. (l45T l) of the standard two-stream instability which is supported by the approximate Eq. (|64] |. Thus 
in Fig. |3]the nearly resonant and the resonant cases are shown with wq — 2ujp and ujq = Aujp, respectively. In the first case the 
high-frequency (w ~ wq — Wp, see Sec. jIVBl l resonant mode is not yet formed but it may interfere with the two-stream mode 
essentially changing the growth rate, see Fig. |3](left panel). And this effect is increasing with the laser intensity. In the second 
case the resonant mode is well separated from the two-stream mode and forms (see Fig. [3] right panel) an isolated maximum at 
^max determined approximately by Eq. ( |59] l (or roughly /cmax — {<^p/c)[t{t — 1)]^/^). The maximal growth rate of the resonant 
mode is well described by Eq. d60] l. It is seen that the position of the maximum is almost independent from the laser intensity 
while the maximum growth rate increases as ~ a^^^ with the parameter a. It is also noteworthy the dependence of the growth 
rate j±{k) of the resonant mode on the frequency loq of the laser field for fixed plasma density and RF intensity. The maximal 
growth rate 7r.max increases with ojq and is shifted towards larger k achieving the maximal value at luq — SAujp. For larger 
frequencies the maximal growth rate is scaled as 7r,max ^ [lLi^p/^o)V ^i^^ f^^^ with ujq. 

Finally in Figs.|4]and|5]the growth rate 7|[ of the longitudinal unstable modes excited in the case Eq || k are shown. These 
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results have been obtained by numerical solution of the dispersion equation (iTZt with Eq. (l73T l for ujq = 1.2ujp, coq = 2ujp 
(Fig.|4]l and ujq = Aujp (Fig.lS). As mentioned in Sec. [V]the coupling between the parametric and the two-stream modes may 
be very effective in this configuration which is supported by the analytical results obtained in Sec. [V] The growth rate (ISSl l of 
the resonant unstable mode in the regime c) is much smaller than 7;^j„ax ^i^^ is not shown here. In both panels of Fig. |4] the 
frequency of the RF is rather small and hence only the right-side resonant mode with 7^,„ax is excited at k+U6 ~ wq + Wp 
(see Eq. (|82]i). The left-side resonant mode with 77,„ax is formed at K_Wb ~ wq — i^p and in Fig. |4]it is merged with the 
two-stream mode and is not distinguishable. To gain more insight in Fig.|5]we demonstrate the growth rates for the larger laser 
frequency luq — Aujp assuming that rib = O.lrie (left panel) and ni, ~ O.Srie (right panel). Now with increasing laser frequency 
Wo the left-side resonant mode is clearly visible in Fig. |5] and the corresponding growth rate 7jrj,,ax is smaller than 7^njax ^s 
predicted by Eq. (|82] |. The domains where the growth rates of the left-side and the right-side resonant modes are nonzero can be 
approximated as Ak± ~ 'i:'yt^max/''^b, see Eq. dSST l. Furthermore, both 7^„jax ^^'^ Ak± increase with electron density as shown 

in Fig.|5](right panel). Note that the growth rates 7^|x ^'^'^ 7r^max increase approximately as ^ n-^ and ^ n-^ with the beam 
density n;,, respectively. In the parameter regimes shown in Figs. 2]and|5]the electron beam is somewhat dense and the condition 
7r^max > 7max i"equii"es rclativistic intensities for the RF (with ve > c) which cannot be fulfilled in the present approximation. 
In the numerical examples shown in Figs. [T]-|5]the parametric instabilities are comparable or even stronger than the standard 
two-stream instability. The stronger effect is expected rather in the case of the transversal polarization Eq J- k of the laser 
field both for the longitudinal and transversal modes. In addition, the fc-domain of the parametric instability is comparable or 
even larger than the range kc of the standard two-stream instability. Although an effective coupling between laser-plasma and 
beam-plasma modes is expected in the case of the parallel polarization of the laser field, Eq || k, the resulting growth rates 
and their /c-domains are in general essentially smaller than those obtained for the transversal polarization. This is not surprising 
because the effects of the laser radiation and the REB on the dispersion properties of plasma are treated here perturbati vely when 
the coupling between both effects is rather weak. 

VII. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have presented a theoretical study of the growth rates of the unstable modes excited simultaneously by 
a laser field and a relativistic beam of charged particles moving in an isotropic plasma. The laser field is treated in the long 
wavelength limit (dipole approximation) and the plasma particles are considered nonrelativistic. In addition, an ultrarelativistic 
beam of the charged particles is considered and the influence of the laser field on the beam is neglected. The dynamics of 
the beam-plasma system in the presence of the RF is studied by the linearized relativistic and nonrelativistic Vlasov kinetic 
equations for the distribution functions of the beam and the plasma, respectively, as well as by the linearized Maxwell equations 
for the electromagnetic fields. The full electromagnetic response of the system is derived in terms of the conductivity tensor 
of the system involving all harmonics of the RF. After a general introduction to the theoretical model in Sec. |II] the dispersion 
relations of the modes are considered in Sees. |III]and|IV] It is shown that in general the longitudinal and transversal modes 
are parametrically coupled due to the presence of the RF. As a result the dispersion equation of the modes represents a secular 
equation for each harmonic of the electromagnetic field. Assuming, however, nonrelativistic laser intensities in Sees. II VI and [Vl 
we have considered only the lowest (zero, first and second) harmonics of the fields truncating the secular dispersion equation on 
the second order of the amplitude of the RF. The dispersion equations derived in these sections led to a detailed presentation, in 
Sec. |VH of a collection of data through figures on the growth rates. For numerical calculations we have chosen ji, — 5 which 
pertains to the FIS relativistic electron beam in the typical 1-2 MeV energy range of practical interest. 

Explicit calculations have been done within cold-fluid approximation both for the beam and the plasma neglecting the low- 
frequency dynamics of the plasma ions. Furthermore, the beam drift velocity is parallel to the wave vector k of the excitations. 
Two particular cases of the transverse (Eq -L k) and the longitudinal (Eq |1 k) polarizations of the RF have been studied in 
detail in Sees. JV] and [V] respectively. For the longitudinal and the transversal unstable modes we have identified some domains 
in the u-k plane corresponding to the resonant and the nonresonant instabilities occurring due to the parametric coupling of 
the different modes. Analytical expressions for the relevant growth rates have been derived which are well supported by our 
numerical calculations in Sec. [VT] These analytical results go beyond those obtained previously in Refs. |43-45] (see also 
Refs. I37ll38ll ) and in Refs. 1 14-16] for the purely parametric and the beam-plasma unstable modes, respectively. In the course 
of this study, we have also derived in Appendix |C] a dynamical equation for the complex amplitude of the excited modes in 
Eo 1 1 k configuration without any restriction on the number of harmonics. When the initial conditions are specified this equation 
may be useful in analyzing the beam-plasma parametric instabilities beyond a weak RF limit considered here. 

It was shown that in the case of the transverse Eq J- k polarization of the RF the longitudinal and the transversal modes 
are coupled due to the RF and can be unstable. The purely beam-plasma and the parametric unstable modes are only weakly 
coupled in this case. With increasing the laser frequency ujq the new branches of the parametrically unstable modes are excited. 
Furthermore, as demonstrated in Figs. [T]-|3]the growth rates of these modes (as well as the corresponding domain in the fc-space) 
essentially increase with luq and the laser intensity and exceed the growth rate 7^^ of the standard two-stream instability. In the 
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case of the parallel polarization (with Eq || k) the longitudinal and the transversal modes are decoupled and the instability occurs 
mainly for the former modes. The purely beam-plasma and the parametric unstable modes are now strongly coupled and the 
whole instability domain in the fc-space is split (at luq > 2ujp) into three major subdomains. In the long wavelength subdomain 
with < fc < fee the instability is similar to the two-stream one while at fcuf, ~ wq ± Wp the instability essentially depends on 
both the beam and the RF parameters (see Figs. |4]and|5). Finally, it was shown that the growth rates are larger in the case of the 
transversal polarization Eq J- k of the laser field. 

Going beyond the present approach, which is based on several approximations and assumptions, we can envisage a number of 
improvements. These include (i) the effects of the low-frequency dynamics of the plasma ions which are completely neglected 
here provided that the obtained frequencies are much higher than the ionic frequencies; (ii) thermal effects both for the beam 
and the plasma. In principle in the case of a standard beam-plasma instability the growth rate can be reduced by these effects 
112811 . However, in the case of a relativistic beam the thermal momentum spreads of the beam particles in the directions parallel 
or transverse to the beam velocity has only little influence on the instability 1 19]. Therefore, these effects are important in the 
case of the nonrelativistic beams as, for instance, in the experiments with heavy ion beams interacting with a laser irradiated 
plasma ll30l432ll : (iii) studying the influence of the finite sizes of the beam in the longitudinal and the transversal directions on 
the instability. An expected effect is the self-consistent generation of the counterstreaming current in a plasma [13] which is 
neglected in the present study. However, this is not a principal restriction on our treatment and the return current can be included 
by adding in Eqs. ( l30b and dSTT i the similar terms but with the flow velocity Ug |13tM1- The latter is determined from the 
condition of the beam current neutralization, rigUe ~ — TifcU;,; (iv) another important issue not considered here is the effect of the 
RF on the dynamics of the beam particles. This implies either relativistic beams (with 7f, 3> 1) or nonrelativistic beams of heavy 
particles (heavy ions, protons, antiprotons etc.); (v) considering the contributions of the higher harmonics with £ ^ 2. From the 
structures of the equations (|39T l. d40] i and dSOt-dTTTi it follows that in this case the combined frequencies iujQ ± ujp come into play 
determining some new branches for the unstable modes. In principle, the study of the generation of the higher harmonics of these 
modes could be facilitated employing a dynamical equation for the complex amplitude similar to Eq. (IC6) ; (vi) lastly, studying 
some other orientations of the polarization vector Eg of the RF and the beam drift velocity u^ with respect to k. However, 
in general, the simultaneous investigation of these issues is a formidable task and requires several separate investigations. We 
intend to address these issues in our forthcoming investigations. 
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Appendix A: Beam-plasma modes in a cold-fluid approximation 



In this Appendix in a cold-fluid model we briefly consider the asymptotic behavior of the frequencies and the growth rates 
of the standard (in the absence of the RF, ve — 0) beam-plasma longitudinal modes at large and small k. Similar qualitative 
analysis have been conducted in Ref. Il28ll . We look for the solutions of the dispersion equation 2?|| (fc, w) = of the longitudinal 
modes in the form uj — uJr + ij, where ujr is the (real) frequency and 7 is the growth rate of the modes, respectively. Here the 
dispersion function Py {k, oj) is determined by Eq. ( |32] i with Eqs. ( |3TT l and (|43] l. 

As has been mentioned in Sec. lIV Al at short wavelengths, k ^ kc, there are four real solutions (i.e. 7(fc) = 0) of the dispersion 
equation which at k ^ kc asymptotically behave as 



UJr2{k) = kUf) 

uJraik) = kub 



1- 



e^^ 



273(fcufc)2 



± 



e^^ 



V?^^ 

'/'few, 2jl/\kubf 



7b KUb 

7^' kub 27^' (kub)^ 



V^u 



+ 



(Al) 
(A2) 
(A3) 



The modes with the frequencies ujr3{k) and w^j^(fc) remain stable also at ^ fc < fee while the other mode with oj^i{k) at 
k — kc merges with the mode u)r2{k). The latter becomes unstable at /c < kc- The frequency and the growth rate of this mode 
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in this long wavelength limit (k < kc) asymptotically behave as 



Wr2(fc) 



kui, 



i + e/7f 






4^(1 -e/7,^) jku.r . 



Similarly one finds the asymptotic behavior of the stable modes i^^i{k) and Wj-a (k) at k < kc, 



:::■«) =T.,yr7^ 



It 



Ui! 



kuh 
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It is seen that the leading terms in Eqs. (IA4l i and (|A5t coincide with LL:J^{k) and 7^^(fc), respectively, see Eq. ( |45T l. 



(A4) 
(A5) 

(A6) 



Appendix B: Evaluation of the sum 

In this Appendix within fluid approximation we briefly derive an analytic expression for the function ^£(C) introduced in 
Sec.|V](see Eq. (iTOll). Inserting Eq. (l43T l into Eq. dTOl l we arrive at 



*£(C) - (-1)' 



E 



S — — 'DC ^ ' 



J,(C)J.+HC), 



(Bl) 



where a± = (w ± Wp)/u}Q. The summation in Eq. ( IB II ) can be easily done using the Newberger's summation formula ll46ll which 
is valid for noninteger /i, ^ 7 ^ 1 and Rc(a + /3) > —1, 



V^ , -,yt '^a-7"(C)'^^+7"(C) _ """ 

^-^ n + /i sin(7r/i) 



•/a+7At(C)>^/3-7A'(C)- 



(B2) 



n— — 00 



In general jjl, a and ^ are complex quantities. Using the summation formula (IB2I > from Eq. (IB II ) at £ ^ and f ^ we obtain 



*KC) = (-1)' 



2wn 



1 



sin(7ra+) 
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(B3) 



2^0 



sin(7ra+) 



Ja^-liOJ-aAO 



sin(7ra-) 



Ja_-^(C)J-a_(C) 



(B4) 



respectively. 

Consider also the limit of the function ^^ {(,) at small parameter C which correspond to the limit of a weak RF. At C ^ 1 using 
the asymptotic behavior of the Bessel function at small argument 114211 from Eqs. ( |B3l l and ( IB4b we obtain 



*^(C) = (-1) 



g+i^i f C, 



w 



E' 

s=0 



(-1)" fe||(-??^s) 

(K|-s)!l + fc||(-7y,s) 



o c 



(B5) 



where r]i — \i\li. 



Appendix C: Dynamical equation for the complex amplitudes 

In Sec. [V] we have derived the dispersion equation for the plasma modes generated simultaneously by the laser radiation 
and the REB. For some applications the derivation of a differential equation for the complex amplitude of the excited wave is 
desirable. We consider here the case of a cold plasma when the partial dielectric function is given by Eq. (|43) . For deriving 
the dynamical equation for the amplitude, we insert Eqs. (|43] i and (|70] i into Eq. (|69] l and using the summation formula (|68] l we 
represent the latter in the form 



ii;ii(o) = - Y. ^iiWfc&iiM) E 



(W - SUJoy - UJ^ 



(CI) 
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Now, in both sides of Eq. ( ICll l we make an inverse Fourier transformation. Then denoting the time-dependent complex amplitude 
as £{t) and using the summation formula I,4Z1 



E 



:Li£uJot 



MO 



±iCsin(tJoi) 



(C2) 



we obtain an equation 

/oo ^ / \2 



(C3) 



Next, on both sides of Eq. ( |C3l l we apply the differential operator d^ /dt^ + ujp. 
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In the last part of Eq. (IC4l i we have used the longitudinal dielectric function of the beam, Eq. (ISTT l. Here Gb{t) is the Green 
function of the electron beam 



Gb{t) 



1 



e'"*da; 



27r y_oo (w — fcufc + i0)2 



(C5) 



where we have introduced the positive infinitesimal +iO which guarantees the causality of the response. An explicit expression 
for the function Gi,{t) can be easily obtained employing in Eq. dCSI ) the contour integration technique. The result reads as 
Gb{t) = 0(— t)te'''"''*, where 8(t) is the Heaviside unit-step function. Substituting this result into Eq. ( IC4b we finally arrive at 



d^ 2 



£ (t) e 



■i{^ sin(cL;oi) 



2 a2 



^b d iCsin(cuot) 
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(C6) 



The obtained equation ( IC6b must be accompanied by the initial conditions. When these conditions are specified Eq 
represents the evolution of the complex amplitude of the waves excited simultaneously by the RF and the REB in a cold plasma 
Note that in contrast to the dispersion equation (|72t . Eq. (|C6l l is also valid for the intensity parameter (^ ~ fca ^ 1. 
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